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ABSTRACT 


An  analytic  performance  evaluation  methodology  for  the  determination  of  the 
expected  value  and  variance  of  the  noise  mean  estimate  of  the  smoothing/prediction 
function  of  the  SALFAS  FM  normalizer  is  presented.  These  statistics  are  necessary  for 
the  determination  of  the  output  Receiver  Operating  Characteristic  (ROC)  curves  of  the 
normalizer. 
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1.0  INTRODUCTION 


The  FM  matched  filter  output  envelope  is  normalized  relative  to  the  local  noise 
background  using  an  estimate  of  the  noise  mean  level  in  each  cell  of  interest  The 
normalizer  algorithm  is  designed  to  provide  robust  and  on-line  adaptive  processing 
features  which  include  [1]: 

a.  Identification  and  removal  from  the  noise  estimation  data  window  interfering 
target  echoes,  unwanted  background  clutter  edges,  and  noise  data  statistical 
outliers  which  would  contaminate  die  window  and  bias  the  noise  mean  estimate. 

b.  Estimation  and  compensation  for  underlying  noise  background  dynamic 
(nonstationary)  intensity  variations  from  cell  to  cell  such  as  may  occur  in 
reverberation-limited  shallow  water  or  convergence  zone  environments. 

c.  A  sufficient  number  of  effectively  independent  samples,  thus  achieving  desired 
random  error  smoothing  and  minimizing  normalizer  SNR  losses  while 
maintaining  desired  control  of  false  alarm  rate. 

Noise  mean  estimation  is  achieved  in  two  steps  as  described  in  [1].  First  the 
matched  filter  output  envelopes  are  operated  on  by  a  sliding  split-window,  median-based, 
peak-shearing  estimation  process.  This  is  designed  to  suppress  interferences  in  the  noise 
data  estimation  window.  The  second  step  uses  a  sliding  split-window  process  with  a 
forward/backward  two  pole  exponential  filter  to  produce  an  estimate  of  noise  mean  of  the 
test  cell  located  at  the  center  of  the  window.  The  data  recursive  filter  allows  for  enhanced 
following  of  arbitrary  noise  mean  level  fluctuations  while  maintaining  desired  random 
error  reduction.  [2] 

Figures  1  and  2  shows  the  median  based  processing  and  the  forward/backward 
smoothing  and  prediction  structures  respectively.  The  window  and  gap  sizes  in  both 
processing  functions  can  be  different  as  shown  in  the  figures.  The  analysis  which  follows 
determines  the  expected  value  and  variance  of  the  predicted  noise  mean  estimates  of  the 
forward/backward  processing  as  defined  in  [1]  which  are  necessary  for  the  generation  of 
the  normalizer  output  Receiver  Operating  Characteristic  curves  [2]. 
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Figure  1  -  Sliding  split  window  structure  for  the  median  based,  two-pass  peak  shearing 
function. 


Figure  2  -  Forward/backward  smoothing  and  prediction  function  window  structure. 
2.0  FORWARD  SMOOTHING  AND  PREDICTION 


The  output  of  the  two  pole  recursive  exponential  smoothing  filter  in  the  forward 
direction  is 


sfh=  rfh-i  +  *  sfh-i] 


(la) 
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#1=  to.,*  a[sfl,  -  tf2 i-i]  (lb) 

where  s/  1,  and  5/ 2,  are  the  outputs  of  filters  1  and  2  at  time  index  i  and  z,  is  the  input 
to  filter  1  (this  is  the  output  of  the  Median  Processing  Function). 

In  Vector-Matrix  form  the  filter  output  is  given  by 

sf.=Asf..  +  bz .  ^ 

—i  —i—i  -  i 

where 

aJ1-*  0 1 

[a(i-a)  i-aj 
£r  =  [a  a2], 

£T  =  Wh  Mil 

and  i  =  {k-G-W  +  1,  ...,k-G-l  }. 


Equation  (2)  represents  the  output  of  the  two  pole  recursive  exponential  filter  at  time  i 
where  i  is  the  sample  index  in  the  forward  smoothing  window  of  size  W,  excluding  the 
leftmost  sample  in  the  window  at  time  k-G-W. 

The  forward  predicted  noise  mean  estimate  of  the  test  cell  at  time  k  in  the  middle 
of  the  normalizer  window  (see  figure  2)  is 

tfk.  -  Zsflk-G-l  '  Mk-G-l  +  — y~  7  (rf h-G-I  -  sf2k-G-l)  ^ 

where  ptfk  is  the  predicted  noise  mean  at  time  k  and  G  is  the  normalizer  half  gap  size. 
The  predicted  noise  mean  estimate  Hfk  is  dependent  on  the  output  of  the  recursive  filters 
at  time  k-G- 1  and  the  constants  a  and  G .  Equation  (3)  can  also  be  written  in  Vector- 
Matrix  form  as 
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(4) 


**  =  £T  U.k-G-1 


.  t  r  t  a(G  +  l)  f,  .  afG  +  l)^ 

where  ci  =  [c;  c2]  and  c,  =  2  +  — — — — ,  c2  =  —  II  +  — — —  L 

Solving  (2)  iteratively  with  an  initial  state  of  sf_kQ  =  lzk0  the  output  of  the  forward 
smoothing  filter  is 

£  =  A‘-"l*k0+  X  A"'  fc  z,  (5) 

tmkO*t 

where  1T  =  [7  7]  and  kO  =  k-G-W. 

The  matrix  form  Am  in  (5)  represents  an  m-fold  product  AxAx...xA(A°  =  I)  and  is 
referred  to  as  the  transition  or  fundamental  matrix  of  a  system.  Inspection  of  equation 
(S)  reveals  that  the  first  term  on  the  right  represents  the  response  of  the  system  to  the 
initial  conditions  only.  The  second  term  represents  the  response  due  to  zero  initial 
conditions  and  the  input  z  and  is  essentially  a  convolution  sum.  Essential  to  the 
determination  of  the  output  is  the  m-fold  matrix  product  Am.  A  brute  force  method  for 
determining  Am  is  undesirable  since  it  would  be  time  consuming  for  large  m  and  insight 
into  the  structure  of  the  filter  response  will  not  be  obtained.  Closed  form  solutions  to  the 
martrix  polynomial  Am  are  obtainable  [3, 4]  and  would  be  shown  later. 

The  expected  value  of  die  forward  smoothing  function  is 

H£i  =  E  {£.}  (6) 

=  A‘-“  !£{*,„)+  £  A-'fe£(z,| 

=  *~k0  l  Htk0  + 

t=kO+l 

where  nt.  is  the  mean  value  of  the  median  processed  data  at  time  sample  i .  Applying 
(4)  the  predicted  mean  in  the  forward  direction  at  time  k  is 
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(7) 


t=kO+l 

W-2 

-  a02k0  +  £ amzk-G-l-m 
m=0 

where  a0  =  (l-a)w  t  (cl  +  c2(l  +  or(IV  -  /))) 

and  am  =  a(l-a)m  (cl  +  ac2(l  +  m)). 

The  expected  value  of  the  forward  predicted  noise  mean  is 

W-2 

£[^]=  a0ntko  +  X amntkClm .  (8) 

m=0 

Utilizing  (7)  the  variance  of  the  forward  smoothing  function  becomes 

aift  =  E(nfl)  -  (») 

=  “o  <>f„  +  2  “o  aT  *  +  aT  R  a 

rw-i  0  0 

where/?  =  0  0  is  a  diagonal  matrix  and 

0  0  rL 

'O’tl-C-,  '/‘Vs-,  rW-2°l  j_c-/  1 

6  _  -G-2°*k-G-l  r°G*k-G-2  rW-3^*k-G-2***k-G-W+l 

•*  “  •  •  •  • 

•  •  •  • 

rW-2^Zk-C-W*l  ®*k-G-l  rW-3®lk-G-W+l  &zk-G-2  r°°*k-C-W+l 

is  a  symmetric  matrix.  The  vector  a  where 

aT  =  [a0  a}  -  otH,_2] 
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is  a  vector  of  the  forward  smoothing/prediction  filter  weights  which  are  a  function  of  the 
filter  constant  a,  the  window  size  W  and  gap  size  G .  The  vector  a  where 

~  **lk_G  2  °*k-G-W+l  ] 

represents  the  standard  deviations  of  the  samples  which  are  mapped  into  the  forward 
window  which  have  been  median  processed.  The  correlation  coefficients  rt_j  are  those 

of  the  median  processed  samples  z,  and  zj . 

The  output  of  the  recursive  forward  smoothing/prediction  filter  (7)  is  the  sum  of 
the  filter  response  due  to  the  initial  state  of  the  filter  and  the  convolution  sum  of  the  filter 
weights  and  the  median  processed  input  The  expected  value  and  variance  of  the  forward 
noise  mean  estimate  are  given  in  (8)  and  (9)  respectively. 


3.0  BACKWARD  SMOOTHING  AND  PREDICTION 

Following  the  same  procedures  as  in  (2.0),  the  output  of  the  backward  two  pole 
recursive  smoothing  function  is 

—j  ~  Asbj+J  +  bzj  (9) 

where  j  =  {k  +  G  +  W  - 1,  ....  k  +  G  +  1)  and  A  and  b  is  as  defined  before.  The 
backward  predicted  noise  mean  estimate  at  time  k  is 


^ k  =  cT  sbk+G+i.  (10) 

Again  solving  for  the  output  of  the  backward  smoothing  filter  iteratively  you  get 

sbj  =  Akl  j  lzkl  +  £  Ai  J  bzt  (11) 

t=kl-l 

where  kl  =  k  +  G  +  W.  The  filter  is  initialized  at  j  =  kl  where  sbkl  =  1 zkI . 
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The  expected  value  of  the  backward  smoothing  function  is 

Us bj  =  E  {skj} 


*Ak‘>  2  £{*„}  +  j  (,£{*,} 


t=kl-l 


=  AkJj  lntkl  + 


t-kl-1 


Using  (10)  the  backward  predicted  noise  mean  estimate  at  time  k  is 

k+G+l 

Hbk  -  cT  Aw  llzkl +  I  A***1  bzt  . 

L  /=*;-/ 


W-2 

=  «0  zkl  +  X  amzk+G+l+m‘ 
m=0 


The  expected  value  of  the  backward  predicted  noise  mean  estimate  is 


TW  — 4. 

4^1-  +  I“»ft,.0.,..- 


The  variance  of  the  backward  predicted  noise  mean  estimate  becomes 


oL  =  E  tnbk2)  ■  E  Ulbt)2 


=  &l  of  +  2  d,  <7,  2r  £  £ 


+  2r  «  2 


where 
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R  = 


rl°t 


*k+G+l 

k+G+2  ^**+G+/ 


^ l^tk+G+l  ^*/k+G+2 

J2 

**+6+2 


r0ai. 


\_rW~2(rZk*G*W-l  °*k+G+l  rW~3a*k+G+W-l  °**+G+2 


rW-2&ik+G*i  ^**+C+W-i 
rH'-30z4+c+2  &tt+G+W-l 

r°G*k+G+W-l 


and  d[ 

zh+G+i+m 


■[« 


**+G+i  **+G+2 


**+G 


•W./]  • 


The  results  for  the  backward  noise  mean  estimate  is  exactly  the  same  as  those  of 
the  forward  noise  mean  estimate  except  that  now  the  samples  arc  those  mapped  in  from 
the  backward  window. 


4.0  FINAL  PREDICTED  OUTPUT  NOISE  MEAN  ESTIMATE  STATISTICS 

The  final  predicted  noise  mean  estimate  at  time  k  is  the  average  of  the  forward  and 
backward  predicted  noise  means  given  by 

Vk  =  |  (tfk  +  Pbk)-  (16) 

The  expected  value  of  the  final  noise  mean  estimate  is 

E  mt]  =  |  (E  lMfk)  +  B  {pbk})  (17) 

where  E{/jfk}  and  E{fibk}  are  given  in  equations  (8)  and  (14)  respectively. 

The  variance  of  the  final  normalizer  noise  mean  estimate  is  given  by  the  sum  of  the 
variances  of  the  forward  and  backward  noise  mean  estimates  and  the  covariance  between 
the  forward  and  backward  noise  mean  estimates. 

The  variance  of  the  final  noise  mean  estimate  is 

<  =  Edit2}  -  Efnk)2  (18) 
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+  2  cov(tfk,tibk)) 


and  the  covariance  between  the  forward  and  backward  predicted  noise  mean  estimates  is 


cov(nfk,nbk)  =  a20alkoalkrkl_k0  +  a0atkl  aT  R  aik  C  hm 

A 

+  &T  *  *  <ZT  *  “ 


(19) 


rW+2G+l  0  0 

where  R  =  0  0  is  a  diagonal  matrix  and 

0  0  r2w+2G-i. 

n 

r2G+2(Tik_G_i  °r*+c+/ 
r2G+3<*zk_G_2  °zk+c+i 

r2G+3  °tk-c-i  °l*+G+2 

r2G+4&zk_c_2  azk+G*2 

•  • 

•  • 

•  • 

r2G+W  c*zk_G_1  CFtk+c+w-i 

r2G+W+l  2  °lk+G+W-l 

• 

• 

r2G+Wazk_c_w+t  <rzk^G+1 

r2G+W +J  ***t+G+i 

r2G+2W-2Ctzk_c_w^J  °*k+G+W-l . 

5.0  MATRIX  POLYNOMIAL  REDUCTION 

The  evaluation  of  the  expected  value  and  variance  of  the  noise  mean  estimate 
requires  the  computation  of  matrix  polynomials  (m-fold  products)  which  can  be  time 
consuming  as  the  order  of  the  polynomial  becomes  large  and  does  not  provide  proper 
insight  to  the  structure  of  the  solution.  Use  of  a  matrix  polynomial  reduction  technique 
can  reduce  greatly  the  number  of  computations  necessary  and  presents  the  matrix 
polynomial  in  closed  form. 


CAYLEY-HAMILTON  TECHNIQUE  [2] 

For  a  matrix  polynomial  N( A)  of  degree  higher  than  A,  if  N(X)  is  divided  by  the 
characteristic  polynomial  of  A  then 
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(20) 


mi 

P(  X) 


=  Q(  X)  + 


m> 

pa) 


where  R(X)  is  the  remainder.  Multiplying  equation  (20)  by  P(X)  you  get 
N(X)  =  Q(X)P(X)  +  R(X). 


If  P(X)  =  0  then  N(X)  =  R(X). 


Let  N(A)  =  Am  so  that  N(X)  =  Xm.  Then  the  characteristic  polynomial  P(X)  of  matrix 
A  is 


P(X)  =  \Xl  -  A\ 


(21) 


=  X2  -  2fiX  +  p2 


where  P  =  1  -  a ,  I  is  the  (2x2)  identity  matrix  and  matrix  A  is  defined  in  (2). 


By  polynomial  division  as  shown  in  equation  (20)  you  get 


\i/ 1 }  _  V* :  oil  im-i-i  .  fnXP  -  (m-l)P 

m>  -vpx  +  +  ' 


(22) 


i=l 


The  first  term  in  equation  (22)  corresponds  to  Q(X)  and  the  second  term  is  the  remainder 
term  corresponding  to  R(X)/P(X).  Since  N(X)  =  R(X)  because  P(X)  =  0,  the 
matrix  polynomial  Am  reduces  to  a  closed  form  given  by 

Am  -  R(A)  =  mpm  l  A  -  (m-1)  pm  I  (23) 


=  (l-a)m 
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6.0  NUMBER  OF  INDEPENDENT  SAMPLES 


The  effective  number  of  independent  samples  in  the  normalizer  window  is  required 
to  be  sufficiently  large  to  achive  the  desired  random  error  smoothing  of  the  normalizer 
estimate.  This  value  is  determined  from  the  ratio  of  the  variance  of  the  normalizer  test 
cell  to  the  variance  of  the  normalizer  noise  mean  estimate.  The  optimum  achievable 
performance  of  a  normalizer  is  its  performance  in  stationary  noise.  Here  performance  is 
determined  solely  by  the  bias  in  the  normalizer  mean  estimate  when  compared  to  the  true 
mean  of  the  test  cell  and  the  variance  of  the  estimate.  The  larger  the  number  of 
independent  samples  in  the  normalizer  window  the  smaller  the  variance  of  the  noise  mean 
estimate  and  therefore  the  better  the  performance. 


The  effective  number  of  independent  samples  is  given  by  the  ratio 


W'/zd 


(24) 


In  stationary  noise  ax.  =  az  and  the  variance  of  the  normalizer  noise  mean  estimate  is 
given  by 


+  rki-ko) 


+ 


W-2 

2<X0  X<*m(rW-7-«  +  rW+2G+l+m)  ’ 

m=0 


q2  W-2  W-2 

j  X  {rm-n  +  r2G+2+m+n )• 

Z  m=0  n=0 


(25) 


7.0  NOISE  MEAN  ESTIMATE  BIAS 

The  bias  in  stationary  noise  can  easily  be  determined.  The  expected  value  of  the 
forward  noise  mean  estimate  defined  in  (8)  is 
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£[#*]=  a0nt  + 


W-2 


m=0 


=  /*« 


«0  + 


W-2 

£a(l-a)“(cl  +  c2(l  +  m)) 


M-0 


} 


(26) 


The  forward  predicted  mean  in  stationary  noise  is  given  by  a  geometric  series.  Using  the 
identities 


W-2 

X/r 


m=0 


1-P 


Y  mfim  =  -  (/  -  /j"'’2  -  (W-2)pw-2  +  (W-2)PWmt) 

£o  (7-/J)2  V  ’ 

equation  (26)  reduces  to 

=  /**  (ci  +  ^2)* 

But  c2  =  1-c,  s  c,  +  c2  =  1.  Therefore  £[#*]  =  Jir.  The  bias  in  the  normalizer 

mean  estimate  in  stationary  noise  is  zero.  Similar  procedure  shows  the  same  is  true  for 
the  backward  prediction. 
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